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Abstract 

The thermal conductivity of the heavy fermion superconductor PrOs4St>i2 was measured down 
to T c /40 throughout the vortex state. At lowest temperatures and for magnetic fields H ~ 0.07ff C 2, 
already 40% of the normal state thermal conductivity is restored. This behaviour (similar to that 
observed in MgB2) is a clear signature of multiband superconductivity in this compound. 

PACS numbers: 71.27.+a, 74.25.Fy, 74.25. Op, 74.45. +c, 74.70.Tx 



The low temperature properties of PrOs4Sbi2, the first Pr-based heavy fermion (HF) 
superconductor |l| (filled skutterudite structure with space group Im3, T c ~1 .85 K) have 
many unusual features both in the normal and in the superconducting states Q]: the non- 
magnetic singlet ground state of the Pr 3+ ion suggests that the conduction electron mass 
renormalization comes from inelastic scattering by crystal field transitions, whereas the 
superconducting transition temperature would be enhanced by the quadrupolar degrees of 
freedom of the rare earth / electrons JjJ |^ . Several experiments also point to unconventional 
superconductivity in this compound: a double superconducting transition in the specific 



heat 0, 0] as well as thermal conductivity measurements in a rotated magnetic field s| 
could result from different symmetry states of the order parameter; London penetration 
depth studies |j| or flux-line lattice distortion [l^] indicate nodes of the gap. Tunneling 



the BCS value, but with a distribution of values 



spectroscopy reveals a gap of the order o 
as observed in borocarbides or in NbSe2 

In this Letter, we report a study of heat transport at very low temperature and under 
magnetic field in PrOs4Sbi 2 , intended to probe the low energy excitations, i.e. the gap 
structure and nodes. Instead, another phenomenon was uncovered: our results provide 
compelling evidence for multiband superconductivity (MBSC) in this compound. Coming 
after similar findings in MgB 2 

H, NbSe 2 Q or in the borocarbides Y and LuNi 2 B 2 C Q, 
they show that very diverse mechanisms may lead to MBSC (or strongly anisotropic gaps) 
so that it could be much more common than presently thought. 

Our rectangular-shaped (~ 0.4 x 0.4 x 2 mm 3 ) PrOs4Sbi 2 single crystal (same as in jsj]) 
was grown by the Sb-flux method |a] and has T c ~ 1.85 K . The thermal conductivity 
(«;) parallel to the magnetic field was measured in a dilution refrigerator by a standard two- 
thermometers-one heater steady-state method down to 50 mK and up to 2.5 T {ixqH c2 {T — > 
0) ~ 2.2 T). The carbon thermometers were thermalized on the sample by gold wires, spot 
welded on the surface of the PrOs4Sbi 2 sample. The same contacts and gold wires were used 
to measure the electric resistivity of the sample by a standard four-point lock-in technique. 

Figure E shows the temperature dependence of k/T at different magnetic fields. Defin- 
ing L = Kp/T, the insert demonstrate the excellent agreement (within 3 %) with the 
Wiedemann- Franz law at the lowest temperatures in the normal state (data in 2.5 T). The 
minimum of L/L at temperatures around 1 K reveals the growth of inelastic collisions on 
warming. For T > 3 K, L/Lq increases above 1, maybe due to a phonon contribution to the 
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heat transport (about 20 % at 6 K). 

At the superconducting transition, k/T (zero field data) exhibits no anomaly, as predicted 
by ordinary BCS theory. The decrease of k(T, H = 0) seems to take place only slightly below 
T c , when the number of excited quasiparticles is reduced by the gap opening. In our case, 
k/T exhibits a significant enhancement at around T c /2 « 1 K. This feature is suppressed 
by a field of only 20 mT (see figure [TJ, whereas the specific heat remains unchanged under 
such small magnetic fields (results not shown here). So this anomaly should be controlled 
by the scattering mechanism. One possibility is an enhanced phonon contribution below 
, suppressed by the mixed state. According to the measured phonon contribution to 
the specific heat and a sound velocity of order 2000 m/s Q], a saturation of the phonon 
mean free path at ~ 20 /im is required to reproduce the temperature and amplitude of the 
maximum of the anomaly. This is 10 times smaller than the crystal smallest length, and 
might come from extended defects in our crystal, like macroscopic voids flux inclusions... 
But other explanations like a boosted quasiparticle inelastic scattering lifetime (as in the 
high-T c cuprates 12, are also possible. In fact, strong enhancement of the microwave 

conductivity is observed below T c , indicating a rapid collapse in quasiparticle scattering 

In the case of a phonon origin of the anomaly, we expect below 0.3 K a contribution to k 
of order 0.2T 3 W/K-m, which should not change the observed temperature dependence of 
our thermal conductivity data below 0.2 K: k ~ 0.26T 2 W/K-m (full line on figure Q). In 
the case of an electronic contribution, it should have disappeared when inelastic scattering is 
suppressed, which is certainly the case below 0.3 K (see data at 2.5 T). The Independence 
of k down to T c /40 indicates low energy quasiparticle excitations. However, it does not 
fit with the simple theoretical predictions for anisotropic gap with nodes: k ~ T 3 for line 
nodes or 2nd order point nodes and k ~ T 5 for linear point nodes. It may result from a 
crossover regime, toward a finite residual value of k/T expected for example in any imperfect 
sample displaying unconventional superconductivity. We also note that experimentally the 
T 2 dependence of k is observed in several unconventional superconductors believed to host 
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231, and CePt 3 Si 



line nodes, such as CeRIn 5 (R=Co,Ir) 21], Sr 2 Ru04 

So let us concentrate on the field dependence k(H) at very low temperatures. As dis- 
cussed above, the quasiparticle mean free path below 0.3 K is governed by elastic impurity 
scattering. We will assume also that the phonon scattering in the same temperature range is 
governed by static defects and is therefore field independent (at least, it cannot be increased 
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by the field). Under small magnetic fields (20 or 100 mT) at very low temperatures (50 
or 100 mK), we observe a pronounced increase in the thermal conductivity with increas- 
ing field (figure H] and |2J) • At intermediate fields, we observe a cross-over to a plateau (for 
H/ H c2 ~ 0.4) which might be related to the symmetry change observed thermal conductivity 
experiments under rotating field 

The if -dependence of k(H) in PrOs4Sbi2 in low fields is in dramatic contrast to that in 
conventional superconductors. For conventional superconductors in the clean limit, small 
magnetic fields hardly affect the very low temperature thermal conductivity. By contrast, in 
unconventional superconductors with nodal structure in the gap function, the Doppler shift 
experienced by the quasiparticles in the mixed state induces a field dependence of k(H). 



The initial decrease of k(H) at high temperature (where the condition yif/if c2 < T/T c is 
satisfied) can be explained by the Doppler shift [25]. But the observed if-dependence for 
PrOs 4 Sbi2 at low fields is intriguing, since it increases with H steeper than that expected 
in the Doppler shift So, though the Doppler shift can explain the low field behavior 

of k(H) qualitatively, it is obvious that it cannot explain the whole if-dependence. The 
extremely strong field dependence of k(H) in PrOs4Sbi2 bears resemblance to that of MgB2 
271 ] (see figure EJ). Indeed, half of the normal state thermal conductivity is restored already 
at H « 0.05if c2 for MgB 2 , and about 40% of k at H « 0.07H c2 in the case of PrOs 4 Sbi 2 . 

As mentioned in the beginning, MgB 2 is now recognized as the archetype of a two-band 
superconductor with full gaps, and it is well established experimentally |i2|, UA l2i| that the 
smallest gap on the minor band is highly field sensitive. Theoretically, for k, as well as for 
the specific heat, the field dependence of the smallest gap is controled by a "virtual" 
(named corresponding to the overlap of the vortex cores of the band with the smallest 
gap (As), having a coherence length of order ^ 0, above H^ 2 , the contribution to « 
of the small band with full gap is close to that in the normal state, only when it is in the 
dirty limit (a condition easily satisfied owing to the large coherence length of that band). 
This remains true even if small inter-band coupling prevents a real suppression of A 9 at 
Hi Q. In the case of P r O S4 S„ 12 , the ,a,e r atio of M m a y originate both f r0m the 
difference in the gap and from the difference in the Fermi velocity between the bands. 

However, this "dirty limit" scenario seems incompatible with the unconventional super- 
conductivity revealed by several experiments ^, U ^| . On the other hand, in case of 
unconventional superconductivity, MBSC might as well give rise to a rapid increase of k(H) 
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at low temperature even in the clean limit. Indeed, for an unconventional superconductor, 
one expects an increase of the contribution of the small gap band on a field scale of order 
Hj? 2 provided the condition ^/H/H^ 2 ^> T/T c is satisfied. In PrOs^Sb^, this will be the case 
at T = 0.05 K or 0.1 K whatever the field above H cl (» 2mT) [3J|. 

An additional experimental observation gives support to the MBSC scenario. Indeed, 
we faced unexpected difficulties in getting reliable measurements in PrOs 4 Sbi 2 compared to 
previous works on other systems (see e.g jsjj): it was not until very thin (17 /jm) Kevlar 
fibers were used for the suspension of the thermometers that reliable values of k (satisfying 
the Wiedemann-Franz law above H C 2) were obtained, and it proved very hard to cool down 
the thermometers below 30 mK. Curiously, if a very small field (~ 10 mT) is applied, 
the thermometers cool down below 15 mK. It is nowadays recognized in the community of 
low temperature thermal measurements, that thermal contacts are a central issue for the 

n 

reliability of such measurements [34| . Before invoking a possible intrinsic mechanism for 
bad thermal contacts (electron phonon-coupling...), we fully characterized these contacts, 
measuring both their electrical (Rf) and the thermal resistance ( R 1 ^ 1 ). 

The results are shown on Figure EJ despite a constant (Ohmic) R e c below T c , R^ di- 
verges strongly at low temperatures, explaining the difficulties encountered on cooling the 
thermometers below 30 mK. It is also seen that this divergence is strongly suppressed in 
a field of 100 mT, much smaller than the upper critical field H c2 , in agreement with the 
observation of the field sensitivity of the base temperature of the thermometers. R e c has 
the usual Maxwell contribution (coming from the concentration of current and field lines 
in the contact area): R% 4 = p/2d, with d — 17 pm the gold wires diameter, and p the 
resistivity of PrOs 4 Sbi 2 |35j (we can neglect the resistivity of the gold wires), which controls 
the temperature dependence of R e c above T c and its jump at T c . It also has an additional 
constant contribution (R cc ps35 mO), coming from scattering at the Au-PrOs4Sbi2 interface. 
We define R 1 ^ 1 as AT/P, P = R e c i 2 the heat power generated by direct joule heating (with 
current i), and AT the thermal gradient across the contact. Following the analysis of R e c , 
AT should have two contributions so that: 

2R% L T/R CC + aT 2 Adn R e c K } 

The first term is coming from AT across the interface. We assume a linear increase of 
the heat power up to R cc i 2 within R cc , a thermal conduction following the Wiedemann- 
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Franz law for the electronic contribution (Lq: Lorentz number), and a aT 2 law for the 
phonon contribution: a ~ 0.18e -6 W/K 3 is the only free parameter of expression and 
has the same value for all fields. The second term is the Maxwell contribution from the 
thermal conductivity of the sample (l/2nd 3^|), with a uniform heat power (R cc i 2 ) plus a 
non uniform heat power generated by R e M (non zero only above T c ) 

Expression gives a fair account of R 1 ^ in 500 mT over the whole temperature range, 
and above 0.8 K in 100 mT and in zero zero field, including the observed jump of R 1 ^ 
at T c (due to the change in the distribution of heat power when R e M is suppressed). The 
additional divergence below 0.8 K might well come from the Sharvin surface resistance, 
which gives additional thermal impedance due to the gap ope ning even in a metallic contact 
(Andreev scattering does not contribute to heat transport, |36|). This divergence of 
whereas R e c remains ohmic and stable put drastic constraints on the thermal insulation of 
the thermometers from the refrigerator, and can be faced in any other experiments with a 
constriction at a normal-superconducting interface. On the other hand, the suppression of 
this divergence in low field in PrOs4Sbi2 is easily explained by the MBSC scenario, because 
above H^ 2 <C H C 2, thermal excitations from the normal metal will be transfered in the small 
gap band without additional barrier. This effect seen on the interface thermal conductivity 
is even stronger than that observed on the bulk thermal conductivity. 

So both k(H) and i?* ft (if), which probe the excitation spectrum, give support to multi- 
band superconductivity in PrOs4Sbi2. A possibility for the origin of multiband band super- 
conductivity in this system is the spread in density of states among the various bands of that 

h A n n 

compound p|: comparison of de Haas- van Alphen and specific heat measurements 
reveals that some of them contain quasiparticles with large effective masses (m* ~ 50m e ) 
and the other only light quasiparticles (m* ~ 4m e ) [stJ, a situation similar to that of most 
Ce heavy-fermion compounds. Theoretical work combining band calculations (for the de- 
termination of vf) and a realistic fit of k{H) (to extract H^ 2 ) is needed to evaluate precisely 
the smallest gap. If instead we take the inflection point at low field of k(H) at 50 mK as a 
"typical value" we get H? 2 ~ 15 mT. With the ratio of the Fermi velocities of both bands 
extracted from H c2 (T) [7j, we find a gap ratio of order 2 : this (rough) estimate contrasts 
with the very large field effect. It is a direct consequence of the hypothesis that the two 
bands of the model have very different renormalized Fermi velocities, which may be taken 
as indicative of weak interband scattering . 
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So our thermal transport measurements under magnetic field provide clear evidence for 
MBSC in the HF compound PrOs 4 Sbi 2 . Strong electronic correlations are thought to be 
at the origin of the different coupling between the various electronic bands, as opposed to 
the difference in the dimensionality of the various sheets of the Fermi surface in MgB 2 . The 
low field behavior of k is consistent with unconventional superconductivity. Further work 
on purer samples with improved thermal contact would be very valuable to precise the low 
temperature behavior of k(T) in this MBSC superconductor in zero field, a study which has 
not been possible in MgB 2 or NbSe 2 owing to a dominant phonon contribution. 
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FIG. 1: k{T)/T at different fields. The anomaly at T ~ 1 K is rapidly suppressed under magnetic 
fields. Full line: pure k/T = aT law valid in zero field below 0.2 K. Arrows: T C {H). Insert: test of 
the Wiedemann-Franz law. Color online. 
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FIG. 2: k(H): Field dependence of k: around T c /2, it may arise from the strong decrease of the 
"IK anomaly" (see figure [J), whereas at low temperatures, the increase signs MBSC (see figure EJ. 
Arrows : H C 2(T). 
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FIG. 3: Low tem per ature behavior of k(H) in a conventional one band 



UPt3 33], MgB2 [22( and our own data on PrOs4Sbi2: adapted from 



superconductors (Nb 



271 ] . The comparison between 



MgE>2 and PrOs4Sbi2 is striking, and supports two band superconductivity in this system. 
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FIG. 4: Squares: R e c of one of the 2 contacts between thermometer and sample (both show the 
same behavior). By contrast, RtJ 1 (circles, here multiplied by LqT) is strongly diverging at low 
temperature in zero field, and highly field dependent. Lines are calculated from expression for 
each field. On this graph, the thermal leak due to the Kevlar suspension remains always above 10 
Ohms. 



12 



